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INTRODUCTION 

Evaluation  of  integrals  involving  generalized  functions  is  often 
accomplished  via  an  integration  by  part-.  w  thout  regard  to  an  interpretation 
of  the  behavior  of  the  function  at  its  singular  points.  See,  for  example, 
ref.  1,  eqs.  (1-17),  (1-18),  and  (1-32);  thus,  there  follows  an  integral  such 
as 


oo 

fdxHiNii  .  _„|W|  .  (1) 

—  *6  X 

It  is  difficult  to  interpret  and  attach  physical  significance  to  this 
integral;  in  fact,  the  major  contribution  to  the  integral  in  (1)  comes  from 
the  neighborhood  of  x  =  0,  where  the  integrand  appears  to  be  positive  and  not 
integrable,  yet  the  right-hand  side  of  (1)  is  negative  and  finite.  We  would 
like  to  approximate  the  generalized  function  1/x2  in  (1)  and  get  a  physical 
interpretation  that  is  consistent  with  the  result  given  by  (1). 

In  a  recent  study  of  the  ideal  patterns  for  arrays  in  one,  two,  and  three 
dimensions,  it  was  found  that  the  required  weightings  were  impulsive  or  more 
singular  than  an  impulse,  depending  on  the  dimensionality  of  the  application; 
see  ref.  2.  In  order  to  make  these  results  of  practical  utility,  it  is 
necessary  to  approximate  tnese  singular  weightings  (generalized  functions)  by 
finite-valued  functions  and  thereby  realize  approximations  to  the  ideal 
patterns.  This  approximation  procedure  and  its  performance  will  be  detailed 
here. 
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APPROXIMATION  PROCEDURE 


We  take  as  given  the  possibility  of  approximating  a  delta  function  £(x) 
by  finite-valued  functions;  see,  for  example,  ref.  1,  page  279,  where  a 
one-siaed  rectangular  pulse  and  a  two-sided  Gaussian  pulse  are  used  for 
illustration  purposes.  Also,  on  page  280,  a  one-sided  approximation  to  the 
doublet  S'(x)  is  given.  Alternative  approximations  are  presented  in  ref.  3, 
pages  11-12. 


Let  us  suppose  that  f(x)  is  an  ordinary  function  which  is  integrable  at 
x  =  a+,  but  that  generalized  function  (ref.  3,  page  30) 

g(x)«f*(x)  (2) 

is  not  integrable  at  x  =  a+.  For  example,  with  a  =  0, 


f(x) 


g(x) 


for  x  >  0 
for  x  <  0 

for  x  >  0 
for  x  <  0 


(3) 

(4) 


is  such  a  pair;  see  ref.  3,  definition  8.  Yet  integration  by  parts  yields 


oo 

J'dx  g(x)  h(x)  =  -jdx  f(*)  h'(x)  =  2  I  dx  x-/4h'(x)  , 


(5) 


which  is  integrable  for  any  good  function  h(x). 


To  make  sense  of  this  situation,  we  approximate  function  f(x)  by  function 
fe(x)  as  shown  in  figure  1;  that  is,  fe(x)  has  constant  value  f(a+e)  in 
the  neighborhood  c  of  a.  The  approximation  to  generalized  function  g(x)  that 
we  adopt  is 
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9c(x)  =  f’(x)  =  f(a+c)  S(x-a)  +  f'(x)  U(x-a-e)  ,  (6) 

as  depicted  in  figure  2,  where  U(x)  is  the  unit  step  function;  that  is, 
9e(x)  has  an  impulse  of  area  f(a+e)  at  x  =  a,  and  is  equal  to  f'(x)  for 
x  >  a+e.  (We  can  approximate  the  impulse  as  discussed  at  the  beginning  of 
this  section.) 

Mx) 

\ 

\ 

\ 

\ 

\ 
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It  is  important  to  observe  that  if  f(x)  is  singular  (ref.  3,  definition 
21)  at  x  =  a,  as  for  example  (3),  the  area  f(a+e)  of  the  impulse  at  x  =  a 
becomes  progressively  1 arger  as  e-*0.  However,  the  area  under  the  remaining 
portion  of  gg(x)  also  becomes  larger,  but  in  such  a  fashion  that  the  total 
area  under  gg(x)  remains  finite  as  e-*0;  this  will  be  shown  below  in  (11). 

To  see  the  effect  of  this  approximation  on  an  integral  of  a  generalized 
function  g(x),  consider  integral 

b 

I  =  f  dx  g ( x )  h(x)  .  (7) 

a 

The  approximation  to  I  is: 


b 

le  s  f  dx  g£(x)  h(x) 
a 


b 

=  f(a+e)  h(a)  +  J  dx  f'(x)  h(x)  ,  (8) 

a+c 

where  we  employed  (6).  Integration  by  parts  on  (8)  yields  two  alternative 
forms: 
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Since  f(x)  is  integrable  at  x  =  a+,  the  leading  term  in  (9A)  and  the  two 
leadinq  terms  in  (98)  both  approach  zero  as  e-*0,  yielding  the  limit 

b 

I  =  1 i m  I£  =  f(b)  h(b)  -  J'  dx  f(x)  h'(x)  .  (10) 

£— ►O  a 

This  is  the  value  of  (7)  expressed  in  terms  of  integrable  functions. 

The  area  under  approximation  g£(x),  for  c  >  0,  is  available  by 
substituting  h(x)  =  1  in  (8)  and  (9): 

b 

J1  dx  ge(x)  =  f(b)  ,  (11) 

a 

which  is  finite  and  independent  of  e.  Thus,  the  impulse  in  figure  2  is 
necessary  in  order  to  compensate  for  the  increasing  area  that  develops  under 
f 1 (x)  near  x  =  a,  when  e  approaches  zero. 

SINGULARITY  AT  x  .  b- 

Suppose,  instead,  that  f ( x )  is  integrable  at  x  =  b-,  but  that  generalized 
function  (2)  is  not.  8y  employing  a  procedure  similar  to  that  above,  we  have 
approximation 

9e (x )  =  -f(b-e)  J(x-b)  +  f’(x)  U(b-e-x)  (12) 

to  g(x).  The  integral  I  can  then  be  expressed  as 
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b 

I£  -  J  dx  ge(x)  h(x) 
a 


b-c 

=  -f(b-c)  h(b)  *  ^  dx  f  ‘  (x )  h(x) 

a 


b-e 


=  -f(b-e)[h(b)  -  h(b-ej]  -  f(a)  h(a)  -  ^  dx  f(x)  h‘(x) 


(13) 


There  follows 


b 


I 


dx  g(x)  h(x)  -  I  =  1 i m  I 

E-*  0 


-f(a)  h(a)  - 


I 


dx  f (x)  h1 (x)  , 


(14) 


in  terms  of  integrable  functions. 

The  area  under  approximation  g£(x)  is  finite  and  independent  of  e: 
let  h(x)  =  1  in  (13)  and  get  directly 

b 

y  dx  ge(x)  =  -f (a)  .  (15) 

a 

RELATION  TO  “FINITE  PART"  OF  INTEGRAL 


If  we  start  with  integral  (7)  and  integrate  by  part  ,  there  follows 


I 


f(x)  h(x) 


b 

^  dx  f(x)  h‘ (x)  . 
a 


(16) 
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Now  if  f(x)  is  singular  at  x  =  a,  the  finite  part  (ref.  3,  page  32)  of  (16)  is 
obtained  by  dropping  the  term  involving  f(a),  thereby  obtaining  identically 
(10).  Conversely,  if  f(x)  is  singular  at  x  =  b,  the  finite  part  of  (16)  is 
just  (14).  Thus,  the  limit  of  the  approximation  procedure  developed  here  is 


exactly  what  is  yielded  by  the  finite  part  procedure. 
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The  result  of  applying  (20)  to  a  good  function  h(x)  is  available  from 
(8)  and  (17), 


X  v+1  p 

1  dx  g£(x)  h(x)  =  h(0)  +  J  dx  x  h(x)  , 


(21) 


or  from  either  form  in  (9).  The  limit  follows  from  (10)  and  (7): 


I  = 


bv+1 

v+1 


r  Vv+1  r 

h(b)  -  l  dx  h'(x)  =  \  dx  x  h(x)  . 

0  0 


(22) 


The  following  examples  are  derived  in  similar  fashion;  just  the  results 
are  listed. 


EXAMPLE  2 


f ( x )  =  In  x  for  x  >  0 


a  =  0 


f ' (x)  =  j  for  x  >  0 


g  (x)  =  In  e  jj(x)  ♦  ~  U(x-c) 


b  o 

I  =  In  b  h(b)  -  J'  dx  In  x  h’(x)  =  J 


dx  i  h(x)  .  (23) 
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EXAMPLE  3 


f(x)  *  j(ln  x)2  for  x  >  0 
a  =  0 


f '  (x)  = 


for  x  >  0 


ge(x)  =  |(ln  c)2  $(x)  +  U(x-c) 

b  b 

I  =  |(ln  b)2  h(b)  -  j”'  dx  ^(1n  x)2  h‘(x)  =  J  dx  h(x)  .  (24) 


EXAMPLE  4 


f(x)  =  Enx-  ^t]  for  x  >  0  ;  v  >  -2 

a  «  0 


f 1 (x)  =  xv  In  x  for  x  >  0 


9,(x)  ■  ['"  *  -  4t]  i(x|  *  x”  x  u(x‘‘> 


v 

^  dx  xv  In  x  i(x)  . 
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EXAMPLE  5 

Here  we  consider  the  result  for  the  generalized  function  1/x2  as 
presented  in  the  integral  (1).  Consider  the  ordinary  function 


f£(x)  s 


-  1/x  for  )x\  >  e 

-  —  sgn(x)  for  |xj  <  e 


(26) 


and  its  derivative 


gc(x)  =  r (x)  =  J 


1/x^  for  jxj  >  e"^ 


(27) 


-  j  £  (x)  for  )x|  <  e 
see  figure  4.  The  function  g£(x)  is  an  approximation  to  generalized  function 
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1/x2;  it  contains  an  impulse  of  area  -2/e  {at  x  *  0)  that  tends  to  as 
e-*0. 

The  approximation  to  integral  (7)  is  as  defined  in  (8);  namely,  for  a  =  -<», 
b  = 


j*  dx  g£(x)  cos(wx)  =  2  J"  dx  C9S(*X)  _  ^ 

<*  6  X 

=  |w|  _  2  1  ~  C-°^I-We?  +  2w  Si  (we),  (28) 

where  Si  is  the  sine  integral  (ref.  4,  eq.  5.2.1).  The  limit  as  e-*0  gives 
the  result  for  the  generalized  function  1/x2,  namely, 


f 


-«0 


dx 


K*  cos(wx) 
yC 


-it  |w|  . 


(29) 

I 


This  result  is  in  agreement  with  ref.  1,  eq.  1-32,  and  with  ref.  3,  page  43, 
for  x-m  with  m  =  2. 

An  alternative  approximation  to  generalized  function  1/x2  that  uses 
finite  functions  is 

exp  (--&)•  1301 


ge(x)  = 


[l lx.6  for  )x|  >  c 
0  for  lx)  <  e 


There  follows 


TR  576; 


f<0 

dx  gg(x)  cos(wx) 


+  2w  Si(we)  .  (31) 

The  limit  as  e-*0  is  again  exactly  (29).  The  essential  feature  of 
approximations  (27)  and  (30)  is  the  large,  sharp  negative  pulse  of  area  -2/e 
that  develops  about  x  *  0  as  e  gets  small. 


=  -it |wj  +  -  cos(we)  -  exp 


(-  f2‘2)] 
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APPLICATION  TO  ARRAY  WEIGHTING 


This  example  relies  heavily  on  material  presented  in  ref.  2  on  array 
weighting  for  good  voltage  response  patterns  in  one,  two,  and  three 
dimensions.  In  particular,  from  ref.  2,  eq.  E-7:  for  -2  <  v  <  -1,  we  have, 

for  the  required  array  weighting  w(x)  to  realize  pattern  v+l^u^“B^), 
one,  two,  or  three  dimensions  0»  or  »  the  relation 


-x  w(x) 


for  0  <  x  , 


(32) 


and  zero  otherwise.  To  match  (2),  we  identify 


g(x)  =  -x  w( x )  , 


(33) 


and 


The  function  f(x)  is  singular  but  integrable  at  x  =  1-,  since  -I  <  v+1  <  0; 
thus  the  derivative  in  (32)  will  generate  a  generalized  function  (ref.  3,  page  30) 
for  the  weighting  w(x). 

We  now  identify  b  =  1  in  (12)  and  obtain,  via  (34)  and  ref.  4,  eq. 

9.6.28,  the  approximation 


9e(x) 


(35) 


Then  there  follows  from  (33),  the  approximate  weighting 
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w  (x) 
€ 


)  S(*-i> 


(¥J  '.(-p)  «,i- 


x)  for  0  <  x  .  (36) 


This  is  the  result  presented  in  ref.  2,  eq.  F— 14.  It  is  an  approximation  to 
the  generalized  function 


w(x)  =  ^2)  Iv(b1T7)u(1-x,  U ( x )  , 


(37) 


which  is  the  required  array  weighting  according  to  ref.  2,  eq.  E-7. 


We  will  carry  this  example  further  than  the  previous  ones,  by  determining 
the  voltage  response  pattern  that  is  actually  realized  by  an  array  employing 
approximation  weighting  (36)  rather  than  (37).  The  array  voltage  response 
pattern  is  given  by  ref.  2,  eqs.  1-12,  for  any  M,  as 

1 

v ( u )  =  J"  dx  x(j)  J^ux)  w( x )  ,  (38) 

0 

where  the  parameter  u  determines  the  array  dimensionality,  and  u  contains  the 
array  geometry,  the  plane-wave  arrival  wavelength,  and  the  various  look  and 
steerir.g  angles.  Substitution  of  generalized  function  (37)  in  (38)  yields 
pattern  (ref.  2,  eq.  E-l) 


v<“>  -Lum 


for  all  u 


where  we  define  entire  function  (ref.  2,  eqs.  17-21) 


(39) 


(40) 
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The  particular  case  of  u+v  =  -1.5  in  (39)  yields  ideal  voltage  pattern 


This  relation  is  true  for  all  u,  whether  larger  or  smaller  than  8.  It 

indicates  a  mainlobe  at  u  =  0  of  amplitude  (§)^cosh(B),  and  sidel^bes  for  u  >  B, 

all  of  equal  amplitude  * 
ref.  2,  figures  13-16. 

When  approximation  w£(x)  in  (36)  is  substituted, instead, in  (38),  we  get 
several  equivalent  expressions  for  the  corresponding  pattern  ve(u),  as  given 
in  appendix  A.  Plots  of  typical  results  are  presented  in  figures  5  and  6  for 
v  =  -1.5.  In  both  figures,  the  curve  labeled  e  =  0  is  the  desired  pattern 
(39).  Figure  5  corresponds  to  a  volumetric-spherical  array  (u  =  .5),  and  the 
desired  pattern  is,  from  (39), 


h\/z  v/  / 

(-1  .  Plots  of  this  ideal  pattern  are  available  in 


which  decays  at  a  3  dB/octave  rate  for  large  u.  Figure  6  corresponds  to  a 
planar-circular  array  (u  =  0)  with  desired  pattern  equal  to  ideal  pattern 
(41).  The  approximations  in  both  of  these  figures  for  e  =  .1  are  quite  good, 
but  those  for  e  =  .2  have  undergone  significant  degradation.  The  possibility 
of  replacing  the  delta  function  in  (36)  by  a  narrow  pulse  is  considered  in  the 
next  section. 


Figure  5  furnishes  an  approximation  to  the  bottom  asterisked  case  in 
ref.  2,  figure  12;  figure  6  does  the  same  for  the  middle  asterisked  case  in 
ref.  2,  figure  13.  The  three  asterisked  cases  for  v  =  -1  in  ref.  2,  figures 
11-13,  merely  require  delta  functions  at  x  =  1  and  are  considered  solved.  The 
last  remaining  asterisked  case  is  the  bottom  one  in  ref.  2,  figure  13,  for 
v  =  -2.  But  this  has  already  been  shown  in  ref.  2,  eqs.  E-37  -  E-39,  to 
involve  a  delta  function  and  its  derivative,  both  of  which  are  easily 
approximated;  see,  for  example,  ref.  1,  pages  279-280,  or  ref.  2,  pages  11-12. 
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USE  OF  NEUTRALIZERS 


An  alternative  smoother  approximation  to  f(x)  and  g(x)  than  afforded  by 
figures  I  and  2  is  by  the  use  of  neutral izers;  see  ref.  5,  section  3.3. 
Consider  the  neutralizer  nc(x)  and  its  derivative  shown  in  figure  7;  the 
neutralizer  is  0  at  x  =  a,  and  1  at  x  =  b.  Furthermore,  it  has  derivatives 


n  (x) 


Figure  7A.  Function  n£(x) 

(x) 


\ 


a +•  b 

Figure  7B.  Function  n^(x) 

Figure  7.  Neutralizer  n£(x)  and  its  Derivative 


of  all  orders,  all  of  which  are  zero  at  end  points  a  and  b.  The  positive 

parameter  e  characterizes  the  critical  point,  xc  =  where  the  neutralizer 
is  1/2;  this  point  xc  will  approach  a  as  e-*Q+.  The  neutralizer  has 
completed  its  transition  to  1  by  the  value  x  =  a-t-e.  n£(x)  will  approach  1 
for  all  x  >  a,  as  e-»0. 

The  approximation  we  take  to  f(x)  is  the  product 

f£(x)  =  f (x)  n£(x)  ,  (43) 

and  the  corresponding  approximation  to  generalized  function  (2)  is  the  smooth 
function 


g£(x)  =  f£'(x)  =  f'(x)  n£(x)  +  f(x)  n£’(x)  .  (44) 

This  approach  is  similar  to  the  regular  sequence  of  good  functions  used  to 
define  a  generalized  function  in  ref.  3,  pages  16-17.  A  representative 
example  is  depicted  in  figure  8.  As  before,  the  area  under  approximation 
ge(x)  is  independent  of  e: 


b 

J  dx  g£(x) 
a 


fe(b)  =  f (b)  . 


(45) 


he  result  of  applying  g£(x)  to  function  h(x)  is  now 


TR  6767 


b 

I£  =  J  dx  g£(x)  h(x) 
a 


i 

J dx  fE'(x)  h ( x ) 

4 


dx  f£(x)  h‘(x) 


b 

=  f(b)  b(b)  -  y 
a 


dx  f(x)  n  (x)  h * ( x )  , 
e 


(46) 


where  we  have  utilized  the  properties  of  the  neutralizer.  Equation  (46)  now 
replaces  (8).  In  the  limit  of  e-*0,  (46)  yields 


U 

I  =  lim  I£  =  f(b)  h(b)  -  J'  dx  f ( x )  h * ( x )  , 


(47) 


e-*0 


in  agreement  with  (10). 


A  limiting  case  of  figure  7  is  a  step-function  at  x  =  an-|;  i.e.. 


n£(x)  =  U(x-a-  |)  . 


(48) 


Then 


nc  * ( x )  =  J(x-a-  |) 


(49) 


and  (44)  yields 


9£(x)  =  f‘(*)  U(x-a-  -I)  +  f(a+-|-)  S(x-a-  ^ ) 


(50) 


This  is  almost  identical  to  (6),  which  has  its  impulse  located  at  x 
rather  than  at  x  *  a*---. 


=  a 
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APPLICATION  TO  ARRAY 


The  starting  point  is  (32),  again,  for  the  required  weighting.  The 
function  f(x)  is  identified  as  in  (34),  and  we  get,  from  (33)  and  (44), 


-xwe(x)  =  ge(x)  =  f*(x)  n£(x)  +  f(x)  n£(x).  (51) 


There  follows,  from  (34)  and  ref.  4,  eq.  9.6.28, 


■  -X  P/ 


f  ■  (x)  -  -x  PV-/  for  0  <  x  <  1  .  (52) 

Substitution  of  (34)  and  (52)  in  (51)  yields,  for  the  approximate  weighting. 


for  0  <  x  <  1 


(53) 


However,  the  neutralizer  in  this  case  must  be  chosen  to  be  1  at  x  =  0, 
and  0  at  x  =  1;  that  is,  it  is  a  reflected  version  of  figure  7A.  For  our 
purposes,  it  is  not  necessary  for  the  neutralizer  to  have  derivatives  of  all 
orders.  Rather,  we  select  n£(x)  so  that  n"(x)  is  continuous  for  all  x, 
and  such  that  n£(x)  and  n£(x)  are  zero  at  the  edges  of  the  transition 
region  (1-e,  1);  see  figure  9.  Here,  letting 


3 


x- 


(r+) 
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Figure  9.  Neutralizer  (55)  and  its  Derivatives 


we  take 


nE(x)  =  j  -  -J(l5y  -  40y3  ♦  48y5)  for  1-e  <  x  <  1 


and  1  for  O^xs  1-e.  Then 


'  i  \  15 

nc(x)  =  -  ^ 


(l-8y2+16y4) 


for  1-e  £  X  i  1, 


and  zero  otherwise.  Also, 


n  (x) 


^ify- 4y3>)  for  1-e  <  x  i  1, 
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and  zero  otherwise.  From  (55)— (57) ,  there  follows 


as  x-*l- 


c  * 


Combining  this  with  the  behavior  of  the  8essel  function  near  zero 
argument  (ref.  4,  eq.  9.6.7),  we  find  that  both  terms  of  the  approximate 
weighting  w£(x)  in  (53)  are  proportional  to  (l-x)3+v  as  x  approaches  1. 

Thus  if  v  >  -2,  the  approximate  weighting  will  approach  zero  at  least  linearly 
at  x  =  1.  An  example  of  w  (x)  for  v  =  -1.5  and  B  =  4  is  given  in  figure  10 


Figure  10.  Approximate  Weighting  w_(x)  for 


for  two  different  selections  of  e.  The  large  positive  spike  near  x  =  1  is 
very  pronounced  for  small  e. 

The  array  voltage  response  corresponding  to  approximate  weighting  wg(x) 
in  (53)  is  given  by  (38)  as 

1 

v£(u)  |  dx  x(i)  J^ux)  w e(x) 

0 

1 

=  J  dx  x2u+1  (ux)  wg(x)  ,  (59) 

0  " 

where  we  used  (40).  A  program  for  the  evaluation  of  (59)  (along  with 
(53)— (56) )  is  presented  in  appendix  B.  Sample  responses  are  plotted  in  figure 
11  for  a  spherical  array  and  in  figure  12  for  a  circular  array.  Comparison 
with  the  corresponding  results  in  figures  5  and  6  reveals  that  the  impulsive 
weightinq  of  (36)  yields  a  better  approximation  to  the  ideal  pattern  (e  =  0) 
than  the  smoother  weighting  of  (53).  That  is,  e  must  be  chosen  smaller  in 
(53)  than  in  (36),  in  order  to  realize  approximately  the  same  voltage  response 
pattern  in  the  first  few  sidelobes. 


Figure  11.  Pattern  of  Spherical  Array  for  Weighting  (53) 
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SUMMARY 

We  have  indicated  how  approximations  to  some  generalized  functions 
generated  by  the  process  of  differentiation  may  be  realized,  and  then  have 
applied  the  procedure  to  the  approximation  of  the  weighting  required  to 
realize  the  ideal  response  patterns  of  arrays  in  several  dimensions.  Two 
examples,  one  impulsive  and  the  other  smooth,  were  used  for  circular  and 
spherical  arrays,  and  the  resultant  approximate  patterns  were  plotted  for 
different  choices  of  the  parameter  e  governing  the  transition  region  near  the 
singularity  of  the  generalized  function.  How  small  e  must  be  chosen  depends 
on  the  form  of  the  approximate  weighting  and  the  desired  closeness  to  the 
specified  pattern. 


29/30 

Reverse  Blank 
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Appendix  A 

VOLTAGE  RESPONSE  PATTERN  FOR  WEIGHTING  (36) 
With  the  shorthand  notation 


approximate  weighting  (36)  takes  the  form 


( A— 1 ) 


we ( x)  =  fv+1(l-e)  J(x-D  +  fv(x)  u(l-e-x)  for  0  <  X  .  (A-2 ) 

Substitution  of  (A-2)  in  (38)  yields,  upon  use  of  (40),  the  corresponding 
voltage  response 


1-e 

V,(u)  =  f^l-e)  >u(u)  +  y  dx  xW^(ux)  fv(x)  ( A— 3 ) 

0 

for  any  u;  this  is  the  relation  programmed  in  ref.  2,  appendix  F. 

For  small  e  (the  case  of  most  interest),  the  cusp  of  f^(x)  in  (A-3)  at 
x  a  1  causes  numerical  integration  problems.  We  can  alleviate  this  problem  by 
integrating  by  parts  on  (A-3),  using 


U  = 


V  =  - 


( A— 4 ) 


to  get 


*.<«>  -  -  <i—)2“jFu(u(i-*>D 


l-e 

+  §  dx  x2u_1  1  _i(ux)  f^jtx)  for  u  >  0  .  (A-5) 


31 


TR  6767 


The  function  fv+j(x)  is  integrable  at  x  =  1  if  v  >  -2. 

The  completed  integral  over  (0,1)  in  (A-5)  can  be  evaluated  in  closed 
form  (ref.  2),  leading  to 


1 

-  J  dx  ■  1#-6) 

1-C 


which  is  advantageous  because  the  interval  (l-e,l)  is  small  for  small  e. 
Finally,  the  cusp  at  x  =  1  can  be  eliminated  for  v  >  -3/2,  by  changing  the 
variable  of  integration  according  to  x  =  sin  t,  getting 


v£(u) 


tr/2 

-  y  dt  (sin  t)2u_1  (cos  t)2'*3  Q  _j(u  sin  <^+1  (8  cos  *  (A_7^ 

A 

where  A  =  arc  sin(l-e),  and 


„  1  U) 

(A_8) 

Form  ( A— 7 )  is  good  numerically  for  small  e  and  any  value  of  u.  There  is  no 
cusp  at  t  *  */ 2  if  \>  *  -3/2.  A  program  for  (A-7)  follows,  where  e,  v,  v,  and 
B  are  arbitrary  (e  >  0,  v  >  -3/2,  B  >  0). 
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18 

Ep*  =  .  1  £  !  Pat  t  ern  '•/*<  u  1  vi  a  ft-7> 

2u 

M  u  =  0  P 

30 

Nu»- 1 . 5 

40 

£  c  =  4  3  f  £  in  1  26 

50 

OUTPUT  O; "Eps  = " ; Eps ;  "  Mu  =  “ ; Mu; "  !U  =  ";Nu: 1 

VM 

II 

60 

DIM  VetO: 480 > 

70 

COM  U, Ec , Ml , N1 , M2, N2 

30 

M  1  =  M  u  -  1 

30 

N 1 =Nu+ 1 

100 

M2=2*Mu- 1 

110 

N2=2*Nu+3 

120 

ft  1  pha=Mu+Nu+ 1 

130 

T=2*Eps-Eps '2 

140 

IF  Eps=0  THEN  160 

150 

F  1  e  =  T  N 1  ■*  F  N I  n  u  x  ri  u  (.  N 1  ,  E  c  *  S  0  R 1  T  >  >  1  ♦"  :•  r  1  '  n  e  2 

«.•  t*  '  ii  -  7  :< 

160 

E 1 = l -Eps 

170 

Te=£  1  ■'  <  2  *Mu  > 

ISO 

A=ASNIE1 > 

1  30 

B=P  I  -'2 

200 

FOR  lu=0  TO  480 

210 

U= I u*P I •  40 

220 

Sq*SQRi AB$<U*U-Bc*Bc ) > 

230 

IF  U  <  =Bc  ThEM  T  1  =FN  I  nuxnu  \  h  1  pt  a,  Sq  !  line 

C  t  .  ri  -  7  ■' 

240 

IF  U)Bc  THEM  T 1 *FN Jnuxnu ( A 1 ph a, Sq > 

250 

IF  Eps : 0  THEN  230 

260 

T3  =  T  1 

270 

V  -  0 

2S0 

GOTO  480 

230 

T2  =  FM Jnuxrtui Mu , U > -T«*FN Jrtuxnu k Mu , U*E 1  ’  1  for  lit 

-.*■  1'  of  V  H ' 

300 

T  3 = T l+Fle*T2 

310 

S»<FNS<fl>+FNS<B) >*. 5 

320 

M  =  2 

330 

H=  <  B-A  )  * .  5 

340 

F  =  <  B-ft> ' 3 

350 

Vo=3E93 

360 

T =0 

370 

FOR  K= 1  TO  M-l  STEP  2 

380 

T=TfFNS<fl+H*K> 

330 

NEXT  K 

400 

S=S  +  T 

410 

V=<.  S  +  T  >*F  !  tins  3  (i  f  v  A — 7  ) 

420 

IF  A  B  S  <  V  -  V  o  >  <  =  A  B  S  <  V  )  *  1  E  -  5  THEN  480 

430 

Vo=V 

440 

N=N*2 

450 

H=H* . 5 

460 

F*F*. 5 

470 

GOTO  360 

4S0 

Ve  <  I  u )“T3-V  !  Vo  1 1  age  Re  sports  e 

h  “  "  > 

430 

PRINT  Iu, Ve< Iu) 

500 

NEXT  Iu 

510 

PLOTTER  IS  "GRAPHICS" 

520 

GRAPHICS 

530 

SCALE  0,480,-70,0 

540 

GRIS  40, 10 

550 

PENUP 

560 

FOR  lu*0  TO  480 

570 

PLOT  I  u ,  20*LGT  <  ABS V* <  I  u  >  ■  Ve  <  0  >  >  ) 

580 

NEXT  Iu 

580 

PENUP 

600 

END 

33 

610 

!  r. 

TR  6767 


1 620  BEF  FHS'T.'  1  i  r>%4gr  an *3  .•»  'H  -  7 

63  0  COM  U,  B: ,  Ml ,  Ml ,  Mi,  N2 

640  St  =  3  I  M  1  T  > 

650  Ct-*COS'T> 

660  T  l=FNJriu.  nu;  Ml  .  U^St 

6 7 0  T  2  =  F  N I r'i u  ' ri u HI  ,  B  c  *  C  t  :■ 

680  RETURN  $t -M2*Ct -'N2*T1*T2 
69  0  F  ME  NIi 

TOO  ! 

710  BEF  F  N  U  aril  m  a 1  X  •  i  ij  i,rn  fn  a  <  X  ■'  '  a  HART,  t:  i  ^  -z  279,  **  3  2'  2*  i 

720  N=  I  NT X  • 

r  SO  R*x-N 

74  0  IF  N :  O >  0 R  '  R  -  ; Cl >  THE N  7 7 0 

750  PRINT  "FNGamma1  IS  NOT  DEFINED  FOR  A  =  * ; 

76  O  STOP 

770  IF  P  0  THEN  SOO 

7  S  0  G  ifn  Hi  a2  =  1 

700  GOTO  040 

300  P«3. 36954359 1  3  1  +  R* I .  09850630453+R*  C .  1 42  92800  71-  49-i-R-*  •  ,  93 

3  1  0  P  *4  3  .  94  1 020’?  1  39+R*  (  22  .  968080082:6 +R  ^ '  12.  SU2  1  6981 1 1  +  P^P  / 

320  0  =  43. 941020'?  191+R*<  4. 39050474596-R*  7 .  1  507503  3299-P  >  > 

320  G*m»*2«P-'Q  1  G»»»»'2+P  ■  for  0  P  <  1 

340  IF  N,-2  THEN  380 

350  IF  N  2  THEN  930 

3  6  O  Gam  m  a=G  am  m  a2 

370  GOTO  930 

830  G  amrna=Gamma2 

390  FOP  K= 1  TO  N-2 

900  G  am  m  a=  G  am  m  a*  X  -  K ) 

910  NEXT  K 

920  GOTO  980 

930  R= 1 

940  FOR  K  =  0  TO  1-N 

950  R=R*<X+K> 

960  NEXT  K 

970  Gamma*  Gam  mat  .'R 

980  RETURN  Gamma 

990  7NEMD 

1000 


0 i 2464 1 
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l  0  1 0 
1  0  2  0 
1  0  5  0 
1  0  4  0 
1 @  5  0 
1  £>(50 
i  o  r  o 
10  20 

I  C  3  0 

I I  oo 
1110 
1  120 
1  1  30 
1  1  40 
1  1  50 
1  1  (50 
1170 
1  1  30 
1130 
1200 
1210 
1220 
1230 
1240 
1250 
1260 

12  70 
1230 
1230 
1  300 
1310 
1220 

13  30 
1  340 
1350 
1360 
1  370 
1330 
1330 
1  400 
1410 
1420 

14  30 
1  440 
1450 
1460 
1470 
1430 
1430 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1530 
1530 
1  600 
1610 
1 620 
1-30 


DEF  F  N  J  r.  u  ri  u  1  N  u ,  !<  >  1  Jr'iU',  \  ri  .i 

IF  hBi  1  THEM  1 1  30 
rl=. 7373S4560SO3 
IF  Mu-0  THEM  RETURN  FNJo-  ■:  ■ 

IF  Mu*.  5  THEM  RETURN  A  +  SIM-'X:-  :: 

IF  Mu*-.  5  THEN  RETURN  H+COS'X' 

IF  N  u  =  1  THEN  RETURN  FN  J 1  ( X  : .  X 
IF  N u  =  -  1  THEN  RETURN  -FN;!  >■ 

IF  Mu*  1.5  THEN  RETURN  A*  <  3 1 N  <  X  > -X*C0S '  X  >  X  3 
IF  H u  =  -  1 . 5  THEN  RETURN  -A*  X*S  I N  < X >  +  C 0 S  ■: X  •  :■ 

IF  N  u  =  2  THEN  RETURN  ■'  £*FN  J  1  X  '>  -X+FN J  o  •  .  3 

IF  Mu* -2  THEN  RETURN  C  2*FN  J  1  X  _■ -X*FN  J  • ..  •  ■  : 

Fi  =  N  u 

IF  '.I  NT  '■  ft  >  A  >  OR  ■  A :  =0  >  THEN  11  63 
k=A=-Nu 

S  =  T  =1  '  2  h»FN 0  a  mm  a.  Ft+  1  >  ;> 

R  =  -  .  254X4* 

El  >3  =  Fl  E  3  '  S  "> 

F 0 P  N *  1  TO  1 0 0 
T  =  T  +  R  '  N  +  ■  N  +  A  •  :> 

S  =  3  +  T 

B  l  g*MAX  Et  g,  AES*'  S >  > 

IF  H  E  3  i:  T  )  *  1 E  —  1  1  *A  £  3  1  3  1  THEN  1  2  7  O 
NEXT  N 

PRINT  "100  TERMS  IN  FN  Jr, u  mu ■  Mu . X  •  AT  •  INuJX 
PAUSE 

D=  1  2-LGT  ■:  HB3  <  E  l  g-'S  )  '•  1  NO.  OF  SIGNIF.  DIGITS 

IF  >  0  THEN  S*3*<4*P>K 
RETURN  3 
FNEND 


D  E  F  F  N  J  o  '■  X  )  !  J  o  <.  X  .>  "i  a  3.4.1  &  3.4.3 

Y  *  A  B  S  <  X  > 

IF  Y  >3  THEM  1330 

T  *  Y  *  Y  -■  3 

Jo*.  0444473-T .  0O33444- T* .  00021  > 

Jo*  1  -  T*  •.  2. 2433337-T*  <  1 . 2656203-T*  ■  . 2 1 6 2366-T- Jo  •  • 

GOTO  1450 
T*3-  7 

Jo =3 . 5 1 2E-5-T+ ( . 00 1 3 7237- T* < . @@@72  S' ©5 -7- . 000 : 44': j  . 
J  o  = .  7  3  7  8  3  4  5  6  -  T  4  (  7’ .  7  E  -  7  +  T  + .  ,  0  U  5  5  2  7  4  0  T  —  J  o  > 

3  = .  0'0262573-T  4  c; .  00054  1  25-r T*  <  .  00023333-7 4 . 000  1  3553  > 
T  =  Y-. 735333 16-T4C. 04 166337+74. r  3. 354E-5-T+S  •  • 

J  o  *  J o 4 C  03  i.  T  )  ■  SQRC  Y  > 

RETURN  Jo 
FNEND 


DEF  FNJl<X>  !  J1<X)  i  a  3.4.4  &  3.4.6 

Y = A  B  S  <  X  3 

IF  Y  >  3  THEN  1550 

T*Y  +  Y.  3 

J 1  =  .  0044  33  1 3-T4 <; .  O002: 1  76  1  -  T  + .  >3000 1  ; 03 

Jl*X+<  .  5-T»<  .  56243335-T4C  .21  03357 3-T  +  .  .  .3 3 35  4  2 3 3 - T  +  J  1 
GOTO  1610 
T  =  3  V 

Jl*. O00 1 7 105-T* .  . 002435 1 1 -T* < . 00 1 1 365 3-T* . 0O020033  1  > 

J  1  *  .  7  3  7  3  3  4  5  6  +  T  4 1 . 5  6  E  -  6  +  T  4  <  .  0 1  6  5  3  6  6  7  +  T  +  J  1  1 

S*  .  03637373-T*C  .  U 0 O 7 4 3 4 3  +  T  +  <  .  000733 24- T * .  @0029;  66  • 

T  =  Y-2 . 356 1 3443  +  T* >: .  1 24996 1 2  +  T4 <;  5 . 65E-5- T  +  S  ,  . 

J 1  =3GN  X  '>  4  J  1  4C0S  < T  ) •-  SQR  <  Y  > 

RETURN  Jl 


FNEND 


"^produced  from' 
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1*40  L  fcl  F  F  N  I  r'i  <j  r*i  m  ■.  Nu ,  .  I  1  1  Irv.,'  r,  '.j 

1*50  if  hB£.':'’  i  them  i  rro 

1**0  ft=.  £ 9 3  94  22 8 0 4 0 1 
1*70  E  =  E  ‘  P  ■ 

1*30  IF  Hu-0  THEM  RETURN  F  M  I  • 

1*90  IF  Nu».5  THEM  RETURN  ft+'E-l  E  ■  X 
1700  IF  Nu=-.5  THEM  RETURN  ft  +  i  E  +  1  E  > 

1710  IF  tliisl  THEN  RETURN  F N I  1  A  ■  X 
1720  IF  Nu=-1  THEN  RETURN  FNIl'I.  +  :’ 

17S0  IF  Nu»1.5  THEN  RE  TURN  ft*<  •  ::-i  X+l  •  £  - .  2 

1740  IF  rk«=-  1 . 5  THEN  RETURN  ft* ■  >  ; I-  1  •  *E-  •.+  !'.£. 

1750  IF  Nu  =  2  THEN  RETURN  0:*FNI*'X  -2  +  FN  I  1  >N  "  £ 

17*0  if  hu*-4  then  return  c ;i»Ftt io -2+fn  i 1  ■ 

1770  ft  =  N  u 

1700  IF  •  I N  T ■  ft  ■  ft'  OR  ■  ft  =0  >  THEN  1800 
1  7  9 0  l=ft=-Nij 

1  *00  *  =  T  =  1  -  ■  2  ft*FNGamma'  ft*  1  >  .* 

1010  P=.25*:>:< 

1820  Big=rtBS<$> 

1*30  FOP  N  =  1  TO  100 
1340  T=T+R.  1 N*<N+ft  ■  > 

1 850  S=S+T 

1  3*0  B  i  g  =  fifty 1  Big,  fit  3  •  S  '  > 

1870  IF  ftB8",  T  >  *lE-li*ftBS'.S  >  THEN  1810 

1880  ne::t  n 

1830  PRINT  "100  TERM  8  IN  FNI  rnj  nuv  rkt,  X  >  ftT  ’  ;  14..4 ; ! : 

1800  PAUSE 

1310  B  =  1  2  -  L  G  T  ■  ft  B  3  ■ Big'S'  >  '  NO.  OF  3IGNIF.  DIGIT: 

1820  IF  K JO  THEN  3=8* < 4*R > - K 
1830  RETURN  3 
1840  FNENB 
1850  1 

18*0  BEF  FNIoOO  1  I o <X  ••  ui  a  3.8.1  s,  3.3.2 

1870  V=ftB8 1  X  > 

1830  IF  Y>3. 75  THEN  2030 
1880  T=Y*Y • 14. 0*25 

2000  I  o  =  .  2  *  5  8  7  3  2  ♦  T  *  •'  .  0  3  *  O  7  *  8  +  T  *  .  0  O  4  5  8  1  3  1 

2010  I  o  =  1  +  T  *  ■  3 . 5  1 5  6  2  2  3  *  T  *  •  3.0899424+T*'  1 . 2  0  *  ‘  4  3  2  ♦  T  *  I  •;> 

2020  GOTO  2070 
2030  T=3.75/Y 

2040  Io=. 008 1*281 -T* ■ . 02057706 -T*- , 0263553"-T*  .01*47*: ?-T*. 00382 377 
2050  I  o=  .  33334228  + T*  <■. ,  01  323592  +  T*  .  00225  318-T*  ■  .  00  1  5 7 5 3  “  -  T  -  I  o  /  >  .» 

20*0  I  c  =  I  o  *  E : ;  P '  .  Y  >  S  Q  R  Y  > 

2070  RETURN  Io 
2030  FNENB 
2080  ! 

2100  BEF  FNI1,:0  i  II •  IO  Mia  3.8.3  4  8.3.4 
2110  Y=ftBS<;  > 

2120  IF  Y  3.75  THEN  2170 
2130  T=Y*Y  14.0625 

2140  I  1  =  .  026537.33  + T*  i  .  00  30 1  532  + T* .  000324  1  1  . 

2150  ll«X*<.5+T*t. 37880594+T *  < . 5 1 498869+T+-  . 1 5034 934 +T- : 1  •  •  ■ 

2160  GOTO  2210 
2170  T  =  3 . 75 - Y 

2130  I  1  =  .  01031555-T*<  .  02282867-T*  •:  .02395312-'  -  .017;  *  :  4-  T  * .  00  »  20059 
2190  I  1  = . 39394228-T * ( . 03933024* T+< . 00 36201 8-T*'  . 00 1  -  2801 - T  *  T i  -  > > 

2200  I  1  =  SGN  CO  *  1 1 +EXP  <  Y  )  /  3QR 1  Y> 

2210  RETURN  II 
2220  FNENB 
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Appendix  B 

PROGRAM  FOR  THE  EVALUATION  OF  (59) 


The  integral  of  interest  is,  from  (59), 


1 

.(u)  =  J  dx  ^(ux) 


x2u+1  W£(x)  . 


( B-l ) 


We  approximate  this  integral  by  sampling  at  increment  a  =  ^  and  using 

"x 

Simpson‘s  rule.  There  follows 


where 


(na)2u+1  w£(na)  , 


( B-2 ) 


{*n}  "  4*  2’  4 . 4*  2’  4’  1  •  <B'3) 

The  program  below  evaluates  (B-2).  Values  of  x2u+^  wc(x)  are  stored  in 
array  Xmuwe. 
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10 

Epi*.  I 

1  TRANSITION  PARAMETER  € 

20 

Hu  =  0 

»  D  J  MENS  I ONAL I  PARAMETER  M 

SO 

Hu*-  1 . 5 

!  WEIGHTING  PARAMETER  y 

40 

Bc=4 

!  WEIGHTING  PARAMETER  B 

50 

N.  =2  -8 

!  NUMBER  OF  INCREMENTS  It)  X 

SO 

N l u*240 

i  NUMBER  OF  INCREMENTS  IN  U 

70  DIM  Xr«u'-.i«  <  0:  1 024  >  ,  Ve<0:  240  ■ 

SO  PEDIM  ::inu»-‘e kO:  N*.» ,  Ve i  O:  Ni  u > 

SO  Mu21=Mu*2+l 

100  No  1 =  Nu+ 1 

110  D* I =  1  N 

120  :i»uuf *. N  >  =0 

ISO  FOP  I=0  TO  H  -1 

i40  x*i*d«i 

150  R2*1-X*X 

ISO  Rs*$QR<R2> 

iro  Br=Bc*Rs 

ISO  I‘'  =  R2'Nu*FNIriUxnu<Nu, Br) 

ISO  IF  X>l-Ep*  THEN  220 

£  O  O  X  m  o  w  e '  I  >  =  X  M  u  2  1  *  I  u 

210  00 TO  2o0 

220  I,,l=R2''Nul*FNInu;<nu<.Nul ,  Br  ■ 

230  N»*FNNs ut < Eps , X > 

240  Nep=FNHeutp<Eps>X> 

250  Xs»uw«<  I  >=X-'Mu21*<  Im*N*-IvJ *N*p-,X> 

260  NEXT  I 

2T0  FOR  I u = O  TO  Ni u 

2 SO  U=Iu/Niu*6*PI 

290  Ud*U*De1 

3  00  T = X  m  u  w  *  <  &  >  +  F  N  J  n  w  x  n  u  > M  u ,  U  *  X  m  u  w  e  ( N :  <  > 

310  So=Se=0 

320  FOR  Ns= 1  TO  Nx-1  STEP  2 

330  So=So+FNJnuxnu<Mu>Ns*Ud> *Xmuw« < Ns  > 

340  NEXT  Ns 

350  FOR  Ns=2  TO  Nx-2  STEP  2 

360  Se*Se+FNJnuxnu<  Mo ,  Ns*Ud )  *Xmuwe  <  Ns  • 

370  NEXT  Ns 

380  Ve<.  Iu)  =  T+4*So+2*Se 

390  PRINT  Iu,V»(Iu> 

400  NEXT  lu 

410  PLOTTER  IS  "GRAPHICS" 

420  GRAPHICS 

430  SCALE  0,Niu,— 76)0 

440  GRID  Niu/6,10 

450  PENUP 

4*0  T  =  20*LGT  <Ve(0) > 

470  FOR  lu=0  TO  Niu 

480  Y=20*LGT < ABS <  W.  I u >  >  > -T 

490  PLOT  Iu, Y 

500  NEXT  Iu 

510  PENUP 

520  END 

530  ! 
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540  DEF  FNNeot  E,:.  j 

556  IF  :><;  =  l-E  THEN  RETURN  1 

560  V>'. 1  ♦ .  5*E  >  E 

570  T  =  Y  *  V 

5SU  Ts  1 5-T* <  40--T*43 '• 

500  RETURN  .5-.125*Y*T 
60U  FNEND 
610  i 

620  DEF  FNNeutp1. E,  X> 

6  50  IF  x:=t-E  THEN  RETURN  0 

64  0  Y  » <  X  -  1  + .  5  •*  E  1  • '  E 

650  T  =  V'4Y 
660  T= 1 -T* < 3- 1 6*T > 

670  RETURN  -1 . 875*T--E 
6S0  FNEND 
6  90  ! 


nLu  OTHER  FUNCTIONS  ARE  LISTED  IN  APPENDIX  fi . 
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